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ABSTRACT 

We present the results of simultaneous XMM-Newton and RXTE observations of the Broad-Line 
Radio Galaxy 3C 111. We find that the Compton reflection bump is extremely weak, however, broad 
residuals are clearly present in the spectrum near the Fe Ka emission line region. When fitted with 
a Gaussian emission line, the feature has an equivalent width of 40-100 eV and full-width at half 
maximum of greater than 20,000 km s~^, however the exact properties of this weak line are highly 
dependent upon the chosen continuum model. The width of the line suggests an origin in the inner 
accretion disk, which is, however, inconsistent with the lack of Compton reflection. We find that much 
of the broad residual emission can be attributed to continuum curvature. The data are consistent 
with a model in which the primary powerlaw continuum is reprocessed by an accretion disk which 
is truncated as small radii. Alternatively, the primary source could be partially covered by a dense 
absorber. The latter model is less attractive than the former because of the small inclination angle 
of the jet of 3C 111 to the line of sight. We consider it likely that the curved continuum of the 
partial covering model is fortuitously similar to the continuum shape of the reprocessing model. In 
both models, the fit is greatly improved by the addition of an unresolved Fe Ka emission line, which 
could arise either in a Compton-thin obscuring torus or dense clouds lying along the line of sight. We 
also find that there are unacceptable residuals at low energies in the MOS data in particular, which 
were modeled as a Gaussian with an energy of ~ 1.5 keV; we attribute these residuals to calibration 
uncertainties of the MOS detectors. 

Subject headings: accretion, accretion disks — galaxies: active — galaxies: nuclei 



1. INTRODUCTION 

In the early 1990s, observations with the Ginga satel- 
lite revealed that the spectra of many Seyfert 1 galax- 
ies contain an Fe Ka emission line, with an equiva- 
lent width (EW) of 100-300 eV, as well as a hard 
excess above 10 keV, relative to the simple power- 
law spectrum fitted over the interval from 2-8 keV 
llPounds et a]|ll990t iPiro. Yamauchi. fc Matsuokalll990t 



IMatsuoka et aLI 119901 INandra Poundijll 994^" These 
features were readily interpreted as signatures of the re- 
processing of the primary X-ray continuum emission by 
nearby Compton-thick material, su ch as the accretion 
disk or the obscuring torus (see, e.g. iLightman fc Whitel 



lips': 'George fc Fabian ^l99li:lMatt. Perola. fc Pirol 
^att et al. 1992). These features are important diagnos- 
tics of the geometry, dynamics, and physical conditions 
of the reprocessing medium. 

In some Se yfert Is, most notably M CG -6-30-15 
ijTanaka et all 119951: llwasawa et all 11996(1 . ASCA ob- 
servations showed that the Fe Kaline profile had a 
narrow core at 6.4 keV and a very broad red wing; 
iFabian et al.l l|1995f) found that the line profile was best 
modeled as emission from the inner regions of the ac- 
cretion disk. In a sample of 18 Seyfert Is observed 
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with ASCA, INandra et al.l l|1997D found that many had 
a broad Fe Ka emission line with an average Gaussian 
energy dispersion of (a) = 0.43 ±0.12 keV. The averaged 
profile of the Fe Ka emission lines, as well as some indi- 
vidual profiles, were not symmetric however, indicating 
that multiple line components (i.e. broad and narrow) 
were present. In some objects the line profile had a broad 
red wing and, like MCG -6-30-15, were well modeled as 
emission from an accretion disk. 

However, recent observations with XMM-Newton and 
have shown that the picture is much more complex. 
Broad lines, typically with equivalent widths of ~ 
100 eV, were detecte d in several Seyfert Is, for exam- 
ple: MCG -6-3 0-15 J Wilms et aLll200l: MCG -5-23- 
16 lIDewangan. Griffith s^fc SchurchI I2003D: NGC 3516 
l|Turner et al."2002 ); Mrk 766 llPounds et alJl200l : and 
Ark 120 fVaughan et all \20m . On the other hand, 
XMM has shown narrow Fe Ka lines (typically unre- 
solved by EPIC) to be a ubiquitous feature of Seyfert 
1 spectra and in fact some Seyfert Is showed onl y a nar- 
row, neutral Fe Kaline, with EW 75 eV I'e.g-. lReevesI 
12003(1 . It must be noted, however, that the upper limits 
for the equivalent width of a broad line were sometimes 
quite generous {EW ~ 100 eV), so a broad component, 
although not required, was not always ruled out by the 
data. 

The hard excess from 10-18 keV observed by Ginga 
is only the low-energy tail of the Compton reflection 
bump, a high energy component which peaks at ~ 30 keV 
and continues up to energies of 100 keV. The strength 
of the Compton reflection bump is parameterized by 
n/2TT, where is interpreted as the solid angle sub- 
tended by the reprocessing material to the primary X- 
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ray source. (In the case of a s t andar d accretion disk, 
f7/27r = 1.) iNandra fc Po^mdsl l|1994D measured f}/27r 
= 0.5-0.7 using Ginga data, but constraining the prop- 
erties of the Com pton bump we l l requ ires very broad 
spectral coverage. iGondek et al.l (^96) combined data 
from Exosat, Ginga, HEAO-1 and GRO/OSSE to ob- 
tain an average 1-500 keV spectrum of 7 Seyfert Is and 
found that f2 /27r = 0.76ib0.15 . More recent observations 
with BeppoSAX (Pero la et aT] l2002: Bianc hi et aL 2001 
and the R ossi X-ray Timing Explorer ( Lee et al]lll999t 
iWeaver ~al. iggS) also indicated that fl/27r ~ 0.7, on 
average. 

Observations with numerous X-ray satellites have 
shown that these reprocessing features are stronger in 
Seyfert Is than in their radio-loud counterparts, the 
Broad-Line Radio Galaxies (BLRGs). The Compton re- 
flection bump i s typically much weaker, w ith n/2TT < 0.5 
(see, e. g.. iZdziarski et all 119951 iWozniak et al.l 
19981 lEracleous. Sa,mbruna. k Mushotzkvl 120001: 
Zdziarski fc Gr andil I2OO1I IGrandi et alT I2001at 
Grandi 2001b). The Fe Kaline is al so weaker, with 
EW 5, 100 eV (see , e.g., ^raclcous. H aloern. &: Livid 
19961 IWozn iak et all 119981 lEracleou s fc Halper nl Il998t 
Eracleous. Sambru n£LfcM ushotz k^i 2000: Grand ieTld] 
2001aHGra ndi 2001b; 'Zdziarski & Grandi 2001). 

Recently, Ballantvnc, Fabian, & Iwasawa (2004) ana- 
lyzed simultaneous XMM and RXTE observations of 
the BLRG 3C 120 and found that Q./2t: ~ 0.5 and the 
Fe Kaline had an EW ~ 50 eV. The line, when fitted 
with a Gaussian, had a width cr ~ 0.1 keV, which is much 
narro wer than those found in Seyfert Is (Nandra ct al. 
119971) . However, a contribution from a broad, distorted 
emission line from the inner accretion disk could not be 
ruled out comple tely. Similar results were obtained by 
iQgle et al.l lj2004(l . using the XMM data only. 

There are several viable scenarios which could explain 
the weakness of the reprocessing features in BLRGs: 

1. The inner accretion di sk in BLRG s might have the 
form of an ion torus ijRees et alJlTgsai . or other 
similar radiativcly inefficient accretion flow. As 
a result, the primary X-ray continuum can only 
be reprocessed by either the outer accretion disk 
or the o bscuring torus, leadin g to < 0.5 
(see, e.g. lEracleous et aLll200nD . In this case, the 
Fe Ka emission line should be narrow (FWHM 
< 15, 000 km s~^) and produced by Fe atoms which 
are not highly ionized (E'--^ 6.4 keV). 

2. iBallantvne. Ross, fc FabianI l|2002fl suggest that 
the reprocessing features are weaker in BLRGs be- 
cause the accretion disk is highly ionized, rather 
than because the geometry of the accretion disk 
is changed. Reprocessing of the X-ray emission 
by ionize d m edia has been studied extensively 
(see, e.g.. iRoss & Fabian 1993; Zvcki ct al. 1994; 
iNavakshin fcKallmanll2001|) . These authors find 
that reprocessing by a moderately ionized accre- 
tion disk results in numerous low-energy emission 
and absorption features, due to ionized species of 
O, C, and N. However, as the ionization increases 
further, the disk becomes a nearly perfect reflec- 
tor, making the reprocessing features very weak. 
In this case, the Fe Kaline should be emitted in 



the inner accretion disk and it should be broad, 
but with > 6.7 keV. 

3. The weak reprocessing features could be the re- 
sult of a mildly relativistic outflow, as sug- 
gested by IWozniak etHl l|1998D . This hypothe- 
sis is supported by detailed modeling of the ef- 
fects of bulk motion on the Fe Ka emission line 
(Rcvnolds & Fabian 199 7jl and the Compton reflec- 
tion bump (;Beloborodovlll99^ . 

In this paper, we present the results of a simultaneous 
observation of the BLRG 3C 111 with XMM and RXTE, 
with the aim of testing the first two of these hypothe- 
ses, which make very clear predictions for the properties 
of the Fe Ka emission line. The general properties of 
3C 111, including the results of previous X-ray studies, 
and the specific goals of this analysis, with respect to 
3C 111, are given In |J3 we describe the data reduc- 

tions. The results of the timing and spectral analysis are 
presented in i^and |Slrespectively. We discuss the impli- 
cations of these results in fj^and summarize our findings 
m in Throughout this paper, we use a WMAP cos- 
mology (Hn = 70k m s"^ Mpc-^^flu = 0.27, f^A = 0.73; 
ISnergel et a].ll200l. 

2. PROPERTIES OF 30 111 

3C 111 is a nearby (z = 0.0485, d = 210 Mpc) 
BLRG. The host galaxy is marginally resolved in the 
R-band with the Hubble Space Telescope and although 
the morphology of the host is somewhat u ncertain, it is 
likely to be a small elliptical-type galaxy ijMartel et alJ 
11999(1 . The ra dio source has an FR II radio mor- 
phologv llFanar off & Rilev 1974) with a single-sided jet 
l|Linfie ld fc Pcrlcv 1984). The jet exhibits superluminal 
motion ( Vermculen fc Cohen 19 94), which along with the 
apparent size of the radio lobes l|Nilsson et al.lll99'^ . al- 
lows us to place constraints upon the jet inclination. As 
described fully in Appendix 1X1 the jet is inclined at an 
angle of 21°-26°. If 3C 111 happens to be a rare Giant 
Radio Galaxy, the inclination angle could be as small as 
10° though. 

A giant molecular cloud lies along the line of 
sight to 3C 111 and some care must be taken to 
estimate the total Galactic Hydrogen column den- 
sity; not only is there a significant contribution from 
molecular Hydrogen, but the molecular Hydrogen col- 
umn density is expected to vary due to the pres- 
ence of AU-scale structures in the molecular cloud. 
Using the H I map of 'Elvis. Lockman, fc Wilke^ 
l)1989l) and detailed H2CO studies of the foreground 
mol ecular cloud by Marsc hCT. Moore, fc Banial l|199^ 
and iMoore fc Marschen (.1995ir the total Galactic col- 
umn density towards 3C 111 is estimated to be 
1.2x10^^ cm~^. This value, is expected to vary by sev- 
eral 10^^ cm~^, however. 

Numerous X-ray satellites have been used 
to obser ve 3C 1 11, mo st importantly Ging a 
iNandra fc Poundsl Il994t IWozniak et alT 119981) 
ASCA (Wozniak et all Il99l IR^eynolds et alJ 11993 
Sam bruna. E raclco ~fc Mushotzkvl I1999D . and RXTE 
(Era cleous et al.l |2()00), which have the high energy 
coverage necessary to study the reprocessing features. 
There are many unanswered questions regarding the 
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X-ray spectral properties of 3C 111 though, that under- 
score the difficulties encountered when analyzing and 
interpreting the weak reprocessing features in many 
BLRGs. 

First, it is uncertain whether a reflection component 
is even necessary to fit the spectrum of 3C 111. Us- 
ing RX TE and Ginqa data resp ectively, lEracleous et al.1 
()2000|) and lWozniak et all l)1998|) found that f2/27r is con- 
sistent with 0, at the 90% confidence level. It is impor- 
tant to place tight constraints upon the continuum emis- 
sion since it has important implications for the geometry 
of the reprocessing medium. Equally important is the 
need to robustly fit the continuum to ensure that the 
residual Fe Ka emission line can be properly fitted. For 
example, when X-ray spectra of 3C 120 were fitte d with 
absorbed powerlaw mo dels (Rcvnolds 1997: Grandi et al.l 
Il997t iSambruna et aD ri999'). the Fe Kaline was found 
to be very strong and broad {EW ^ 0.5 - 1 keV 
and cr 1 keV). However, when Com pton reflection 
or a broken powerlaw models were vised ( Wozniak et al. 
[1998; Eracleous et al. 2000; Zdziarski fc Gr andi 2001), 
the line width and the equivalent were significantly re- 
duced {EW - 100 eV and a ~ 0.3 keV.) 

Secondly, while the equivalent width of the 
Fe Ka emission line in 3C 111 was found to be weak 
{EW ~ 60 tw eV) b y both Eracleou s et al. i2000) and 
I Wozniak et al.l lll998D. the energy and w idth of the line 
were uncertain. lEracleous et all l|2000ri fitted the line 
with a Gaussian with a fixed energy of 6.4 keV and con- 
strained the full-width half-max ( FWHM) of the line to 
be less than 44,000 km s"\ while iWozniak et alJ l|1998l^ 
fitted the line with a narrow Gaussian {a = 0.1 keV) 
and found that E = 6.7 ^q's keV. Furthermore, if the 
Fe Kaline originates in the accretion disk, a Gaussian 
model is a poor approximation to the true disk line 
profile and yields misleading estimates of the line 
energy and width. The li ne was margina l ly det ected 
in an ASCA observation ()r!,evnolds et al.l I1998D and 
the line properties were not well constrained. Thus, 
neither the origin of the line (i.e. inner disk vs. outer 
disk, or an even more distant reprocessor, such as the 
torus) nor the ionization state of the reprocessor are 
known. In order to evaluate the competing scenarios 
to explain the weakness of the reprocessing features 
of BLRGs presented in ^ these parameters must be 
better constrained. 

A simultaneous XMM-Newton and RXTE observa- 
tion of 3C 111 can help address these and other is- 
sues. It is important to obtain simultaneous obser- 
vations, since the X-ray flux and spectral parame- 
ters of AGNs in general, and BLRGs in particular, 
are known to vary on timescales of several days (e.g., 
|gI iozzi. Sambruna fc Eracleou'^ l2003j) . The high en- 
ergy sensitivity of RXTE is critical for accurately fit- 
ting the continuum and detecting the Compton reflec- 
tion bump, which peaks at 30 keV. On the other hand, 
the good spectral resolution and large collecting area 
of XMM-Newton in the 0.4-10 keV range make it ideal 
for use in a detailed study of the Fe Ka emission line 
properties. Additionally, the spectrum at low ener- 
gies will be useful for constraining the ionization of 
the disk because reprocessing by a moderately ionized 
disk leads to numerous emission and absorption features 



at low energies llRoss fc FabianllT99l IZvcki et al.llT99l 
iNavaks hin & Kallmanll2001fl which should be detectable 
in the XMM-Newton spectrum, despite the large absorb- 
ing column. 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. XMM-Newton 

3C 111 was observed on 14 March 2001 with the Eu- 
ropean Photon Imaging Camera (EPIC) and the Reflec- 
tion Grating Spectrometer (RGS) on-board the XMM- 
Netwon satellite for a duration ~ 40 ks. The p-n data 
were obtained in Large Window mode and the MOS data 
were taken in the Partial Window mode, using the thin 
filter. The exact exposure times and count rates for each 
instrument are given in Tabled The data were processed 
using the XMM-Netwon Science Analysis Software (SAS 
v5.4.1) using the calibration files released on 29 January 
2003. The EPIC data sets (p-n, MOS 1 and MOS 2) 
were filtered to remove all flagged events (e.g. events sus- 
pected to be cosmic rays, bad pixels, etc.). There were 
intense particle flares during most of the second half of 
the observation in which the count rate increased by a 
factor of 20-100. The flares were successfully removed us- 
ing the Good Time Interval tables provided, which also 
removed several smaller flares that were present. As a 
result, the exposure times were reduced by ~ 40% (see 
Tabled]). 3C 111 is quite bright and normally the effect of 
the flares might have been adequately corrected for with 
background subtraction. However, we noticed that the 
source count rate actually decreased dramatically during 
the flaring intervals, leading us to believe that the flares 
must have been intense enough to saturate the telemetry, 
despite the fact that the count rates were well below the 
expected threshold for telemetry saturation (the satura- 
tion thresholds are 1150 and 300 counts for the p-n 
and MOS detectors respectively.) Therefore, the data 
collected during the flaring intervals could not be used 
in any way. Finally, the EPIC p-n data were filtered to 
include only single and double pixel events (i.e. PAT- 
TERN < 4) whereas the MOS data were filtered to also 
include triple pixel events (i.e. PATTERN < 12). 

The source counts were extracted from a circle with a 
radius of 44" for the p-n data, for which the fractional en- 
circled energy is greater than 90%. The background was 
extracted from an annulus centered on the source with 
an outer radius of 110" and inner radius of 45"; we exper- 
imented with several background regions and found little 
difference in the results. Since the MOS data were ob- 
tained in Small Window mode, the largest circular region 
which could be used to extract the source spectrum had 
a radius of 39", which encompasses 88% of the encircled 
energy. The background could not be extracted from the 
same chip, because it was not read out. Instead a circle 
of radius 125" was extracted from a neighboring chip. 
As with the p-n data, several measurements of the back- 
ground were found to be similar. The response matrices 
(RMFs) and ancillary response functions (ARFs) were 
generated with the calibration files released on 29 Jan 
2003. When making the ARFs, the source was treated 
as a point source, but the encircled energy was modeled 
as a function of photon energy. 

The RGS data were reduced using the SAS routine 
RGSPROC, which automatically filtered the data, traced 
the 1st and 2nd order dispersed image of the source, and 
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selected regions to exclude in the determination of the 
background spectrum. The RGS data suffered from the 
same flares as the EPIC data, however we found that 
it was unnecessary to remove the flares; the background 
subtracted spectra with and without the flares removed 
were similar. When the spectra were extracted, a sepa- 
rate background file was created for each order. 

In the preliminary spectral analysis, it became clear 
that the p-n data suffered from X-ray loading, which oc- 
curred because the frames used to calculate the offset 
map (i.e. the zero-energy level) had a count rate which 
was too large. The offset map, and thus the gain, were 
therefore incorrect on a pixel-by-pixel basis. Addition- 
ally, because the zero-energy level was too high, double- 
pixel events were interpreted as single-pixel events of a 
higher energy. In many cases X-ray loading is an extreme 
effect of photon pile-up. However, we note that based on 
the count rate of the filtered p-n data, the chance of pho- 
ton pile-up is < 1% both from the estimates from the 
XMM-Newton handbook and our own estimates based 
on Poisson statistics. The X-ray loading in this instance 
likely was the result of a flare in one or more of the frames 
used to calculate the offset map. At this time, there is 
no method to reliably correct for X-ray loading. Thus 
we were forced to exclude the p-n data from the spectral 
analysis. With the total loss of the p-n data and the 40% 
loss of MOS data, due to flaring, the total count rate was 
reduced by 70% from that anticipated. 

The X-ray loading manifested itself in several ways, 
but was difficult to diagnose. There were two revealing 
symptoms of X-ray loading that we noticed, which distin- 
guished the phenomenon from pile-up. First, when the 
p-n and MOS spectra were fitted with simple absorbed 
power-law models, there was an absorption feature be- 
tween 1.8 - 2.1 keV in the p-n data residuals which was 
absent in the MOS data. We initially suspected a calibra- 
tion error in the effective area near the Si edge, however 
this feature was absent in other p-n data with similarly 
high signal-to- noise ratio (S/N). Secondly, the results of 
the SAS EPATPLOT routine indicated that the fraction 
of single pixel events was slightly higher than expected 
from theory whereas the fraction of double pixel events 
was slightly lower, the opposite of what was expected 
for photon pile-up. However, this effect was not dra- 
matic and might have been easily overlooked had there 
not been suspicious residuals in the p-n data near 2 keV. 

3.2. Rossi X-ray Timing Explorer 

3C 111 was observed with the Proportional Counter 
Array (PCA) and the High Energy X-ray Timing Exper- 
iment (HEXTE) onboard the Rossi X-ray Timing Ex- 
plorer (RXTE) satellite from 14-17 March 2001. The ex- 
posure time for each instrument is given in Table [T| The 
reduct ion procedure is described in detail in lGliozzi et al.l 
l|2003|) . Briefly, the PCA and HEXTE data were screened 
to exclude events taken when the Earth elevation angle 
was < 10° and the pointing offset from the optical po- 
sition was > 0.02°. The PCA data were also filtered to 
include only events obtained when the satellite was out 
of the South Atlantic Anomaly for more than 30 min- 
utes and also those events whose ELECTON-0 param- 
eter was < 0.1. The PCA background and light curve 
were determined with the L7-240 background developed 
at the RXTE Guest Observer Facility (GOF), using the 



FTOOLS task pcabackset, v2.1b. The appropriate re- 
sponse matrices and effective area curves for the obser- 
vation epoch were produced using the FTOOLS v. 5.1 
software package and with the help of the REX script 
provided by the RXTE GOF. Only PCUs and 2 were 
combined, since PCUs 1,3, 4 were not always turned on. 
The background applicable to the HEXTE clusters was 
obtained during the observation by dithering the instru- 
ment slowly on and off the source. 

4. TIMING ANALYSIS 

To perform the timing analysis of the XMM data, the 
MOS 1 and MOS 2 data are combined and the p-n data 
are used, as a slight offset in the gain should not affect 
the timing results. The data from 0.2-10 kcV are binned 
in 2000 s intervals. The mean count rate of the p-n data 
is 8.89 s"^ and that of MOS 1-^2 is 6.98 s"^ The count 
rate is moderately variable, with an amplitude of 2.7%, 
where the amplitude is defined as the difference between 
the maximum and minimum count rate, divided by the 
mean count rate. However, the variability is significant 
and the probability that the count rate is constant is 
only 2.1%. The data are consistent with a monotonic 
increase in flux during the course of the observation, but 
the hardness ratio (defined as the ratio of the flux in 
the 2-10 keV band to that in the 0.2-2 keV band) is 
not significantly variable, and the probability that it is 
constant is 43%. The 2-20 keV data from the PCA are 
also binned in 2000 s intervals and the mean count rate 
is 17.29 s~^ for 2 PCUs. The variability is consistent 
with the XMM data in the overlapping interval with an 
amplitude of 2.5%. As can be seen in Fig. over the 
entire observation period, the lightcurve is more variable 
and the probability that the count rate is constant is 
< 10"'^*. The largest excursion has an amplitude of 11% 
and takes place over a 29 ks time interval. As with the 
XMM data, the hardness ratio (defined as the ratio of 
the flux in the 10-20 keV band to that in the 2-10 keV 
band) is not highly variable and the probability that the 
hardness ratio is constant is 82%. 

5. SPECTRAL ANALYSIS 

In total we have nine separate data sets, obtained 
simultaneously, covering the spectral range from 0.4- 
100 keV: RGS 1 and 2 (0.4-1.65 keV in first order 
and 0.65-1.65 keV in 2nd order); MOS 1 and 2 (0.4- 
10.0 keV); PCA (4.0-30.0 keV); and HEXTE clusters 
and 1 (20.0-100.0 keV). The RGS data are binned 
such that each bin has a minimum of 25 counts. The 
MOS data have an excellent S/N, therefore they are 
binned such that there are 2-3 bins per resolution ele- 
ment, and at least 25 counts per bin. In particular, the 
region around the expected Fe Kaline is well sampled. 
At low energies (E < 0.7 keV), the MOS resolution ele- 
ment was slightly undersampled, because the count rate 
is not large. However, this region overlaps with the high 
resolution RGS data. The RXTE data were binned such 
that there were at least 20 counts per spectral bin to en- 
sure that the test was valid. The spectral analysis is 
carried out with the XSPEC vs 11.3 software package 
(Arnaud 1996). 

The data were obtained simultaneously, so the fit pa- 
rameters for the nine data sets are forced to be the same 
in all models, with the exception of the overall normaliza- 



RXTE and XMM observation of 3C 111 



5 



tion constants which are allowed to vary freely to account 
for cross calibration uncertainties. In general, there is 
only a 1% difference between the normalization constants 
for the two MOS data sets, but the PCA normalization 
constant is 30% higher than the MOS value. Since the 
RXTE observations span a longer time period than the 
XMM observations, it is possible that spectral variability 
will lead to a systematic difference in model parameters 
between the XMM and RXTE data, a possibility which 
we explore below. 

All errors and upper limits listed in the Tables and the 
text correspond to the 90% confidence interval for 1 inter- 
esting degree of freedom (d.o.f., i.e. Ax^ = 2.7), unless 
otherwise stated. In comparing different models, we refer 
to the chance probability Pc, by which we mean the prob- 
ability that the improvement in the fit statistic would 
occur by chance, as determined by the F-test. In ii5.ll 
we describe fits to the continuum using various models, 
excluding the energy interval where the Fe Ka line is ex- 
pected. The residual Fe Ka line from several continuum 
fits is modeled in H5.2I Then in ^b.'3\ we attempt to fit the 
combined continuum and Fe Ka emission simultaneously 
and self-consistently. 

5.1. Continuum Models 

We fit the continuum using several different mod- 
els, excluding the interval from 4.5-7.5 keV, where an 
Fe Ka line is expected to be located. The fit parameters 
for all models are listed in Table El and discussed below. 
In all models, we include Gala ctic ph otoclectic absorp- 
tion u sing the cross-sections of iMorrison fc McCammonI 
l) 198.3ft . allowing the column density to be a free param- 
eter because the total Galactic Hydrogen column den- 
sity along the line of sight is uncertai n. Throughout, we 
use th e solar abundance pattern of lAnders fc Ebiharal 
()1982i) and in all instances the fitted column density is 
~ 8 X 10^^ cm~^, which is consistent with the Galactic 
Hydrogen column density towards 3C 111.^ Therefore 
an additional absorber at the redshift of 3C 111 is not 
warranted, although absorption within the host galaxy 
or the AGN itself cannot be ruled out. 

Powerlaw (Model #1): - The data are first fitted with 
an absorbed powerlaw model, whose free parame- 
ters are the column density [Nu ) , the photon index 
(r), and the overall normalization constant. The 
spectral model and residuals are shown in Fig. El 
and the 90% confidence contour in Nh and F is 
shown in Fig. 

We use this simple model to test the assumption that 
there are no systematic differences between the XMM 
and RXTE data sets. Because F and Nh are only loosely 
constrained by the RGS and HEXTE data, we do not 
include those data sets in this test. We allow Nh and 
F for MOS 1 and MOS 2 and F for the PCA data to 
vary independently, but we set A^//for the PCA data 
equal to that of MOS 1 because TV// is unconstrained 
by the PCA data alone. The MOS 1 and 2 data sets 
yield consistent values of Nh , but there are discrepan- 
cies in F, with Fmosi = 1-74 ± 0.02, Fmos2 = l-78t°o*l?2 

^ A similarly good fit is o btained using the abundance pattern 
of lAnders fc Grevessd il98HI) , but only if the Oxygen abundance is 
reduced by ~ 50%. In this case, Njj ~ 10^^ cm^^. 



and FpcA = 1-69 .q q^. We find that this discrepancy 
is not the result of the spectrum, as whole, becoming 
harder during the longer RXTE observation; fits to sev- 
eral temporal subsets of the PCA data are best-fit with 
the same value of F. When the MOS 1, MOS 2, and 
PCA data are fitted above 3 keV (excluding the interval 
between 4.5-7.5 keV), we find that Fmosi = 1.57 
FmoS2 = 1.68 ±0.08 and FpcA = 1.70 j;|^:^3. Thus it ap- 
pears that the discrepancy in the powerlaw index is an 
indication that the spectrum hardens at high energies. 
However, it is clear that the photon indices inferred from 
the MOS 1 and PCA data are inconsistent, although the 
MOS 2 and PCA data are in excellent agreement. This 
can be seen in Fig. (21 the MOS 1 data show clear positive 
residuals above the Fe Ka line region that are not present 
in the MOS 2 data. Throughout the rest of this paper, 
we will assume that the XMM and RXTE data can be 
fitted with the same set of parameter values and allow 
only the normalization constants to vary independently, 
however the effects of allowing the MOS 1 data to have 
an independent photon index will also be considered. 

Broken Powerlaw (Model #2): - The next simplest 
model is an absorbed broken powerlaw, whose free 
parameters are the column density, two powerlaw 
indices (Fi and F2), the breaking energy (Eh), and 
the normalization constant. The fit is greatly im- 
proved, with a chance probability of ~ 10~^^. The 
small breaking energy, E^ — 2.9 ± 0.2, suggests 
that the addition of a low energy component might 
improve the fit greatly. However, this extra com- 
ponent also increases the complexity of the model; 
before including it, we first investigate the effects 
of including Compton reflection from neutral and 
ionized media. 

Powerlaw + Compton Reflection from a neutral medium 
(Model #3): - We next fit the data with a model 
which includes powerlaw emission from the primary 
X-ray source as well as the continuum from X-rays 
reprocessed by a cold, dense slab of material. This 
model, i mplemented with the PEXRAV routine in 
XSPEC ijMagdziarz fc Zdziarskil IT9951) . does not 
include Fe Ka line emission nor does it include the 
gravitational effects which would arise if the repro- 
cessing material is near the black hole. The free pa- 
rameters are: Nh , F, the folding energy of the pri- 
mary X-ray powerlaw continuum (Enfold), the cosine 
of the inclination angle (cos i < 0.95 or i > 18°), 
and the solid angle subtended by the disk to the 
prima ry X-ray source {n/2Tr). Like[Eraclcous et aQ 
we find that the elemental abundances do 
not greatly affect the fit and we keep them fixed to 
the solar values. As described in Appendix fXl the 
inclination angle can be further constrained to be 
less than 26°. As E'foid increases above 850 keV, 
the improvement in the fit is negligible, so we re- 
strict E'foid < 1000 keV. With these additional con- 
straints in place, we find that there is a marginal 
improvement in the fit as compared to the simple 
powerlaw (Pc = 17%). The spectral model and 
residuals are shown in Fig. [3 and the 90% contour 
in Nh ~F space is included in Fig. As shown 
by the Eioid-^/^n confidence contours, (Fig. 0, 
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ri/27r is small, but non-zero. We note that f2/27r 
is insensitive to the inclination angle, within the 
narrow range of allowed values. 

Powerlaw + Compton reflection from an ionized 
medium (Model #4): - As discussed in the 
ionization state of the disk is important in de- 
termining the strength of the reflection features, 
so we now consider reflection from an ionized 
disk, as implemented with the pexriv rou tine in 
XSPEC (Magdziarz & Zdziarski 1995; Don e et all 
I1992D . The additional free parameters in this 
model are the disk temperature, T (< 10^ K) 
and the ionization parameter, ^ = 47r Fion/n, 
< 5000erg cm s~^) where Fion is the incident 
5-20 keV flux and n is the density of the reflect- 
ing medium. The improvement in the fit, as com- 
pared to a neutral reflection model, is negligible 
(Pc = 37.6%) and we found that < 20erg cm s~^ 
and n/2T: = 0.13 ±0.06. 

When the above models are fitted to the data, there 
are unacceptable residuals at low energies (See Fig. 
Furthermore, the results of the broken powerlaw flt in- 
dicate that the addition of a soft component is likely to 
improve the continuum flt significantly. Thus, we now 
add a variety of low energy components to a powerlaw 
continuum in an attempt to improve the flt. 

Powerlaw + Thermal Bremsstrahlung (Model #5): - 
We flrst add a thermal Bremsstrahlung component 
to the absorbed powerlaw model, leading to a sig- 
nificant improvement in the flt compared to the 
powerlaw model (Ax^ = 67 for three extra d.o.f.; 
Pc ^ 10^^"^). The thermal component has a tem- 
perature of ~ 1 keV and is found to contribute 
10% and 1% of the flux in the 0.5-2 keV and 2- 
10 keV bands, respectively. The spectral model 
and residuals are shown in Fig. |21 and the 90% 
contour in Nh -F space are included in Fig. |31 We 
note that a ~ 0.3 keV blackbody component also 
provides a good flt. Even with the addition of the 
Bremsstrahlung or blackbody component however, 
some low-energy residuals remain. 

Powerlaw + Soft Gaussian (Model #6): - It is also pos- 
sible that the low-energy residuals are the result of 
remaining problems in the calibration of the MOS 
effective area curves. Therefore we also attempt to 
model the residuals with combinations of emission 
lines and absorption lines and edges. We find that 
a single Gaussian emission line with ^ 1.5 keV, 
a ~ 0.5 keV, and a flux of ^ 10^^^ erg cm^^ s^^ is 
quite effective in eliminating the majority of the 
low energy residuals (Ax^ = 89 for three extra 
d.o.f; Pc ~ 10~^^). The spectral model and resid- 
uals and the 90% contour in Nb -F are shown in 
Figs. 12 and 13 respectively. 

Whether a ~ 1 keV Bremsstrahlung, a ~ 0.3 keV black- 
body, or a Gaussian emission line is added, the best-fit 
powerlaw photon index decreases to the value we find 
from fitting only the PGA data (F = 1.69). Therefore, it 
is not surprising that when these components are added 
to the models which include reprocessed emission, no 



Gompton reflection is necessary and il/2TT < 0.05. We 
find that if we exclude the data from 1-3 keV, where 
the low-energy residuals are the most extreme, we find 
again that F = 1.69; it appears that the addition of a low 
energy component is truly necessary to obtain a reliable 
description of the 3-100 keV spectrum. The origin of 
this component is discussed further in ii6.1l 

In summary, the data are best fitted with a simple 
Powerlaw (with Nh ^ S x 10^^ cm"^ and F « 1.7) with 
little or no contribution from Gompton reflection of the 
primary powerlaw emission. This is in marked contrast 
with Seyfert galaxies, in which {Q,/2tt) = 0.7 (see fQ. 
Throughout the remainder of this paper, a soft Gaussian 
component is included in all models because it provides 
the best statistical fit, however the results we obtain us- 
ing a thermal Bremsstrahlung (or blackbody) component 
are similar. 

5.2. Models for the Fe Ka Line 

In this section, we study the profile of the residual 
Fe Ka emission from the Powerlaw -t- Soft Gaussian and 
Powerlaw -I- Gompton Reflection models (#6 and #3 re- 
spectively). Before adding the line component, all data 
outside the 3.5 - 8.5 keV energy interval are excluded and 
the continuum parameters, including the normalization 
constants, are also frozen. The emission line is fitted by 
either a Gaussian line or a relativistically broadened line 
from an acc retion disk, implem ented with the diskline 
routine (see iFab lan et al.lll989fl . The Gaussian parame- 
ters are the line energy (E), the energy dispersion (cr), 
and the line flux. The diskline model has the follow- 
ing parameters: the inner and outer radius of the line- 
emitting region (Rm and Rout, where R is expressed in 
units of the gravitational radius, rg — GM,/c^, with Af, 
the mass of the black hole), the disk inclination (z), the 
slope (3 of the emissivity pattern, (e oc r^), the energy of 
the emission line (E), and the line flux. Note that the 
diskline model does not include the effects of the Gomp- 
tonization of the Fe Ka line as it emerges from the disk, 
which asymmetrically broadens it further. The line pa- 
rameters for all three data sets are forced to be the same, 
with the exception of the normalization constants, which 
are linked together such that the equivalent width of the 
line was the same for all three. 

Since the various line model parameters are highly in- 
terdependent, we show the flt results in the form of con- 
tours (Figs. ISlEl andlTJ, rather than list them in a table. 
For the Gaussian model, 90% confidence contours in the 
line energy and the FWHM of the line profile and con- 
tours in the EW and the FWHM are shown. For the disk 
line model, 90% confldence contours in the inner radius 
of the disk (i?m) and the line energy are constructed for 
models with /3 = -2, -2.5 or -3 and i = 10, 18 or 26°, 
with Rout flxed to be 500 rg ; the contours are very similar 
for Rout = 1000 rg. ^ 

The improvement in the fit with the addition of an 
emission line, whether modeled by a Gaussian or a disk 
line, is significant. For the Powerlaw -I- soft Gaussian 
continuum model — 314 for 172 d.o.f. when no line is 
included, whereas = 163 for 169 d.o.f. for the Gaus- 
sian model. For the disk line model, = 159, 165, and 
173 for 169 d.o.f. for /? = -3, -2.5, and -2 respectively 
and Ax^ is not very sensitive to the inclination angle. 
For the Powerlaw -I- Reflection continuum continuum. 
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— 279 for 172 d.o.f. when no line is included, whereas 

— 184 for 169 d.o.f. for the Gaussian model. For the 
disk line model, = 186, 187, and 190 for 169 d.o.f. 
for (3 = ~3, -2.5, and -2 respectively. The residuals are 
equally well modeled by the Gaussian and the disk line, 
but for the Power law + soft Gaussian continuum model, 
a disk line model with /? = -2 is disfavored by the data 
because of the large value of x^- 

As discussed in i i5.1l and seen in Fig. 13 the MOS 1 
data show clear residuals above the Fe Ka line region 
because those data prefer a smaller F than the MOS 2 
and PCA data. Because the properties of a weak line 
are very dependent upon the continuum fit (see §2), we 
also present the results for a fit in which a Gaussian 
line and a simple powerlaw continuum were fitted si- 
multaneously only over the region from 3-10 keV. The 
photon indices were allowed to vary independently but 
Nh ^B.s fixed to be 8 x 10^-^ cm~^, because it is uncon- 
strained by these high energy data. The parameters of 
the Gaussian line were initially allowed to vary indepen- 
dently. We found that the line parameters inferred from 
the MOS 1 (i^MOSi = 6.3 ± 0.2 keV, ctmosi < 0.5 keV, 
and EWmosi = 50 ± 40 eV) and PCA [EpcA = 6.4 1°;^ 
keV, crpcA < 0.5 keV, and EWpca = 50 tto eV) data 
were in agreement. The line is not well constrained by the 
MOS 2 data, but when E and a are the same as for MOS 
1, the upper limit to EW is 70 eV, which is consistent 
with the EW measured for the MOS 1 and PCA data. 
Thus for the purposes of investigating the line properties, 
the line energy, a and EW were forced to be the same 
for all three data sets.^ As with the Powerlaw -I- soft 
Gaussian and Powerlaw -I- Reflection continuum models, 
the addition of a narrow line leads to a significant im- 
provement of the fit (x^ = 155 for 220 d.o.f. compared 
to 210 for 221 d.o.f when the line is not included.) 

An Fe Ka emission line is clearly required to fit the 
data regardless of the continuum model. The line is re- 
solved in all cases with a large bulk velocity (FWHM 
> 20, 000 km s~^ and Rin < lOOr^ for the Gaussian and 
disk line models respectively) which implies an origin in 
an accretion disk as opposed to a distant reprocessor such 
as the obscuring torus. However the line need not form in 
the innermost regions of the accretion disk. The equiva- 
lent width also has considerable uncertainty {EW = 40- 
110 eV), but is much smaller than the total EW of the 
narrow -I- broad Fe Ka emission line (EW ~ 200 eV) ob- 
served in those Seyferts with broad Fe Ka emission lines. 

5.3. Combined Continuum and Fe Ka Emission Models 

In the previous two sections, we verified that reprocess- 
ing features, especially the Compton reflection bump, 
are much weaker in 3C 111 than in Seyferts. How- 
ever, the models used are not physically reasonable; the 
Fe Ka emission line is very broad, implying an origin 
in the inner accretion disk, but the Compton reflection, 
which would necessarily accompany this disk line, is not 
allowed by the high energy continuum. In the follow- 
ing sections, we simultaneously fit the continuum and 

^ Although the line is not resolved by any individual data set, 
when all three are fitted together the line must be broad, as seen 
in Fig. 5. This is because the high S/N PCA data require that the 
line has EW > 40 eV, but for an unresolved line the upper limit 
to the MOS 2 EW is 30 eV; to fit all three data sets, the line must 
be broad. 



Fe Ka emission line, with the goals of finding a self- 
consistent model for the data and assessing whether ei- 
ther of the two scenarios presented in ^ are physically 
reasonable. The soft Gaussian component is included in 
all models and the best-fit parameters are listed in Table 
Eland discussed below. 

5.3.1. Truncated Accretion Disk - Models #7a,b 

The total spectrum from 0.4-100 keV is fitted with 
a model which includes Compton reflection and an 
Fe Ka emission line, both of which arise from the re- 
processing of the primary X-ray emission by a neutral 
or moderately ionized accretion disk (Model #7a). We 
use the refsch routine, which implements the ionized 
disk model described in ^5.1\ (Model #4) , but convolved 
with a disk line model to ac count for relativistic effects 
(see, for example iFabian et al. 1989). The Fe Kaline 
is modeled as a disk emission line with /3 = —3, i = 
26°, Rout = 500 rg and the normalization of the emission 
line is tied to fi/27r such that EWpcKa = 160 x fl/27r eV 
(jGeorge fc F abian 1991), which is only appropriate for 
a solar abundance of Fe. Both the continuum and line 
parameters are allowed to vary. 

The inner radius of the accretion disk is Rin ~ 
llOl^Q^Tg and the energy of the emission line is con- 
sistent with being emitted by nearly-neutral Fe, thus 
this model indicates that both the continuum and 
Fe Ka emission line could arise from a truncated accre- 
tion disk. However, there are residuals in the spectrum 
near 6.4 keV. Usi ng the Chandra High Energy Trans- 
mission Gratinef. lYaaoob fc P admanabhan (2004) found 
that the narrow cores of Fe Ka emission lines in Seyferts 
have an average energy of ~ 6.4 keV and a ~ 0.02 keV. 
The inclusion of such a narrow Gaussian, in which E and 
cr are fixed to these values, improves the fit slightly (Ax^ 
= 8 for one d.o.f.. Model #7b). The equivalent width of 
this n arrow line is ~ 25 eV. Ghiselli ni. Haardt. fc Mafr3 
(jl994|) found that for disk inclinations similar to that of 
3C 111, a Compton-thin torus (with Nh ^ 10^^ cni^^) 
will contribute an Fe Ka emission line with an equivalent 
width of 30 eV while making a negligible contribution to 
the Compton reflection bump. The line can also origi- 
nate by transmission through dense clouds along the line 
of sight. However, with the addition of the narrow line, 
the reflection fraction, and consequently the equivalent 
width of the emission line arising from the disk, are very 
smaU (r2/27r = Omta°Ql). Given the weakness of the 
disk line it is impossible, with these data, to place any 
meaningful constraints on either the energy of the disk 
emission line or the inner radius of the accretion disk. 

The results of these flts suggest that only a very weak 
{EW ~ 10-20 eV) broad emission line is required and 
the data are primarily fitted by a narrow emission line. 
This is in sharp contrast to the results of ii5.2l which 
suggest that the Fe Kaline is stronger {EW — 40-110 
eV) and very broad {FWHM > 20, 000km s~^). This 
is due in part to the fact that there is significant broad 
curvature in the refsch continuum. Using the fakeit 
command in XSPEC a data set was simulated based 
upon the refsch model parameters for Model #7a, but 
without the disk emission line component; when these 
data are fitted with the Powerlaw -I- soft Gaussian model 
(#6), there are broad residuals in the Fe Ka emission line 
region, as seen in Fig. |S1 This curvature is likely due 
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to the relativistically blurred Fe Kaedge, as we noticed 
that the excess emission was greater for larger values of 
J7/27r but less noticeable for smaller values of Rm- If 
the DISKLINE component of Model #7a is replaced with 
a Gaussian, the line only has EW = 40 ± 10 eV and 
a = 0.3 ^q'^ keV, which is not inconsistent with the total 
EW (narrow + broad) of ^ 40 eV inferred from model 
#7b. Thus the broad residuals studied in ii5.2l can be 
fitted with a combination of continuum curvature as well 
as broad and narrow line emission. 

The fit results for this model (#7b) are consistent with 
the emission arising from the reprocessing of the pri- 
mary X-ray continuum by an accretion disk, although 
a distant reprocessor also contributes to the flux of the 
Fe Kck emission line. Because the component of the line 
arising from the accretion disk is very weak, Ri„ is not 
constrained and the data do not require the disk to be 
truncated. However a larger value of Rin (such as found 
for Model #7a) is more physically attractive because it 
is consistent with the small value of n/2Tr found. 

5.3.2. Highly Ionized Accretion Disk - Models #8a,b,c 

As discussed in ^ the reflection features from a 
highly ionized disk are very weak; when fitted with 
the PEXRiv and refsch routines in XSPEC, as we 
have done thus far, the value of n/2TT is mislead- 
ingly low ijRoss. Fabian, fc Youn e 1999). Therefore, 
we now fit the data with the constant density ion- 
ized disk model described by lRoss fc Fabianl (| 1993( ) and 
iBaUantvne. Iwasawa. fc Fab ian (2001), which is avail- 
able as a table model for use in XSPEC (Model #8a). 
We refer to this model as RF-PEXRIV to avoid con- 
fusion with the previously discussed pexriv model. 
IBaUantvne et alJ (|J002ll compared this model to the 
ASCA spectrum of the BLRG 3C 120 and found that 
the spectrum is well fitted with ^ = 4000erg cm s~^ and 
r2/27r fixed at unity. 

The disk is modeled as a slab of gas with a constant 
density of 10^^ cm~^, which is illuminated by a powerlaw 
with index F and a sharp cut-off at 100 keV, and the ion- 
ization parameter extends beyond 10^. (Here the ionizing 
flux is defined from 0.01-100 keV.) This model includes 
the emission from the Fe Ka line, so the emission line and 
continuum are fit self-consistently. The Comptonization 
of the reprocessed photons as they emerge from the dense 
disk is also included, but gravitational and inclination- 
dependent radiative transfer effects are not. 

Even though the best-fit ionization is quite large (^ ~ 
4000) the reflection fraction is Q/2Tr = 0.3 ±0.1; n/2TT ~ 
1, as observed in Seyferts, is not allowed (see Fig. EJ. 
As before, the addition of an narrow emission line with 
E = 6A keV improves the fit slightly (Ax^ = 7 for 1 
extra d.o.f; Model #8b); the EW of the narrow line is 
~ 15 eV but the reficction fraction is further reduced to 
n/2TT = 0.2. The results of these fits suggest that even 
when the disk is allowed to be highly ionized, ^l/2^T is 
still rather small. However, if gravitational effects are 
included (as in §5.3.2), larger values of n/2Tr might be 
allowed, because the refiection features would be blurred. 

As a test, we have convolved the rf-pexriv model 
with a disk emission line^ with /3 = —3, i = 26°, Rout = 

® The convolution was carried out using a code kindly provided 
by A.C. Fabian. 



400 Tg and also included a narrow (unconvolved) 6.4 keV 
Gaussian (Model #8c). The response functions were ex- 
tended using the extend command in XSPEC and the 
HEXTE data above 80 keV were ignored because the 
model is only defined up to 100 keV. We find that the 
fit is improved when the relativistic effects are included 
but that the best-fit parameters (see Table 3) are very 
similar; in particular fJ/27r ~ 0.2 and EWn ~ 15 eV. 
(A model in which i?i„ = lOOr^ fits equally well.) This 
model is computationally expensive and it is impossi- 
ble with the available computer facilities to extensively 
explore parameter space for this model or to determine 
error bars for the parameters. Instead, we determined 
best-fit models for specific values of r2/27r (0.6 and 1.0) 
and we find that these larger values are not absolutely 
ruled out, but they are not favored by the data (Ax^ — 
+7 and -1-13 for 1 d.o.f. respectively, compared to the 
best fit model #8c). 

The models presented in this section provide very good 
fits to the data and it is certainly possible that the accre- 
tion disk is ionized. However, the second scenario pro- 
posed in ^ (that U/2tt ^ 1 and the reprocessing features 
are weak only because the disk is highly ionized), is not 
preferred by the data, although it is not completely ruled 
out. 

5.3.3. Partial Covering - Model #9 

Recent observations with XMM have shown that, con- 
trary to the early results from ASCA, many Seyfert 
Is possess only a narrow Fe Ka line with no evidence 
for a broad, distorted line (Reeves 2003). In some of 
these objects, broad residuals are clearly present in the 
Fe Ka region of the spectrum, but they are equally well 
fitted by a broad diskline and a model in which the con- 
tinuum is partially covered by an absorber with a column 
density of lO^^-ss ^,^-2 (ggg^ g g_ 1H0419-577, Pounds 
et al. 2004a; NGC 4051, Pounds et al. 2004b). This 
heavily absorbed powerlaw spectrum turns over at an 
energy near the Fe Ka emission line, thereby introducing 
curvature into the continuum spectrum than can mimic a 
broad Fe Ka line. Here we test whether the broad resid- 
uals in the Fe Ka region of the spectrum of 3C 111 are 
simply an artifact of partial covering. 

The data are fitted with an absorbed powerlaw model 
that includes an additional partial covering absorber with 
column density iVH,2 and covering factor /c, implemented 
with the PCFABS routine. There are clear residuals in 
the spectrum so a Gaussian emission line is also included 
(Model #9). Confidence contours for iVH,2 and fc are 
shown in Fig. 1101 The emission line has an energy of 

6.4 keV and an EW ^ 30 eV and, as can be seen in 
Fig. ^2 the line is not resolved at the 90% confidence 
level. We find that even when the data from 3-8 keV are 
excluded from the fit, the same values of iVH,2 and Fc 
are found, which suggests that continuum itself is bet- 
ter fitted by this partial covering model than a simple 
powerlaw. We note that this model, in which none of 
the primary emission is reprocessed by a dense, geomet- 
rically thin accretion disk, could be considered to be an 
extreme case of scenario #1 presented in fll (i-^- the 
accretion disk is truncated). 

6. DISCUSSION 

6.1. Origin of the Low Energy Component 
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Many of the continuum models presented in HS.ll do not 
provide a satisfactory fit to the low energy data. We find 
that the inclusion of a low energy component, described 
either by a ^ 1 keV Bremsstrahlung (or ~ 0.3 keV black- 
body) or a broad Gaussian line with E ^ l.h kcV and 
(T 0.5 keV, greatly improves the fit. Although it ap- 
pears that this component must be included to obtain a 
good description of the 3-100 keV spectrum, it would be 
reassuring to determine whether there are any plausible 
sources of this low-energy component. 

Previous ASCA and i?05'^T observations have shown 
no evidence for a soft component in the spectrum of 
3C 111, however the large and uncertain column density 
would have made its detection difficult. Soft components 
have been observed in other BLRGs. The BLRG 3C 382 
shows a soft excess which cannot be entirely explained by 
the known e xtended halo of hot ga s which surrounds the 
host galaxy l|Grandi et al.ll2001a|l . The BLRG 3C 120 
exhibits a clear soft excess, that contribute s 20% of the 
0.6-2 keV flux that 'Ballant vne et alJ l|200# model with 
a ~ 0.3-0.4 keV Bremsstrahlung. However, the origin of 
these soft components are not well understood. 

An alternative possibility is that the soft energy com- 
ponent simply compensates for an error in the calibra- 
tion of the MOS effective area curves. The count rate 
near 1 keV is very large, with some data bins contain- 
ing more than 4000 counts. The statistical uncertainties 
(~ 2%) are on the order of the uncertainties in the cal- 
ibration of the MOS effective area curves, which can be 
as large as ~ 10% (XMM-Newton helpdesk, private com- 
munication). For this reason, we also modeled the low 
energy component as a Gaussian with i? 1.5 keV and 
a ^ 0.5 keV. The peak amplitude of the emission line 
is only ~ 8% of the continuum flux, so it is not impos- 
sible that the soft Gaussian component is compensating 
for an error in the MOS effective area curves. The soft 
Gaussian component is admittedly ad-hoc. However, this 
component has been chosen, instead of a thermal com- 
ponent, primarily because it affects a relatively isolated 
region of the spectrum (compared to the Bremsstrahlung 
model) while at the same time it produces a better fit 
to the low energy residuals. There is independent ev- 
idence that there is a known excess in the MOS data, 
similar to what ar e present in these data (see Fig. 11 in 
iKirsch et al.ll200'^ . We have verified that including the 
Soft Gaussian component only in the two MOS data sets 
does not affect the results presented here. 

6.2. Interpretation of the Spectral Models 

As seen in ii5.1l and ii5.3l the high energy data do not 
allow for a significant Compton reflection bump, even if 
the disk is highly ionized (see i)5.3.2l Models #8a,b,c). 
Therefore these data disfavor the presence of a standard, 
optically-thick, geometrically-thin accretion disk (neu- 
tral or ionized), which extends to small radii, in which 
case Q,/2tt ^ 1. For comparison, if the outer, thin ac- 
cretion disk transitions at some radius to a vertically ex- 
tended structure, (such as a n ion-torus), then < 0.3 

l|Ghen fc HalT)ernl[l 989: Zdz iarski et all 19991) . consistent 
with the values of Sl/27r we have obtained. 

Although the results of ii5.2l suggest that there is a 
broad {FWHM > 20, 000 km s^^ Fe Kaline in 3C 111, 
these residuals can be partly attributed to curvature in 
the continuum, either due to relativistic blurring of the 



Fe Kaedge f ii5.3.1(l or partial covering by a dense ab- 
sorber fi i5.3.3|l . Weak Fe Kaline emission with EW ~ 40 
eV is still required, but the majority of the flux in this 
line most likely arises in a distant reprocessor. Thus, 
from the point of view of fltting the spectrum of 3C 111, 
it is not necessary to invoke reprocessing by an accretion 
disk to model the data and there is no obvious reason to 
prefer the truncated disk models (#7a,b) over the more 
simple partial covering model (#9). 

However, it is not immediately obvious which struc- 
ture could be obscuring the primary X-ray source. The 
inclination of 3C 111 is less than 26° so unless the ob- 
scuring torus has a very large opening angle or is very 
extended vertically, it cannot be responsible for the par- 
tial covering. An alternative source of partial covering is 
Compton-thin {Nh < 10^'^ cm~^) clouds along the line 
of sight. These same clouds could also contribute to the 
unresolved Fe Ka emission line. The expected EW of the 
Fe Kaline transmitted by such clouds is (Halpcrn 198^ 

350^ (6A\ / N,., \ /_4^\ 
T + 2\7.lJ yiO^^cm-^ J \4 X 10-y 

where fc, Nh,2, and F are the same covering factor, 
column density, and powerlaw index as found from 
the partial covering model flt, and Ap^^ is the Fe 
mass fraction. The above equation reduces to EW 
= 71/c(iVjj/1023 cm~2) gv. For a given EW of the 
Fe Ka emission line, there is one-to-one relationship be- 
tween Nh,2 and fc] tracks for EW = 10, 15, and 20 
eV are shown in Fig. ^| As can be seen, the par- 
tial covering absorber can contribute an Fe Ka line with 
an equivalent width of at most 15 eV. However, if the 
clouds have a super-solar abundance, the Fe Ka emission 
would be enhanced. It is likely that the bulk of the 
narrow Fe Ka emission line is formed in the obscur- 
ing torus, which can contribute ~ 30 eV to the to- 
tal equivalent width. Transmission through clouds with 
Nh,2 > 10^^ cm~^ could also make a small contribu- 
tion to the observed Fe Kaline, but it is very difficult, 
without detailed modeling, to estimate the EW of an 
Fe Ka emission line in this case. 

Although BLRGs, on average, do not have large intrin- 
sic column densities ijSambruna et al.|ll999() . the BLRG 
3C 445, which also has an FR H radio morphology, is 
similar to 3C 111 in several respects. The intrinsic ab- 
sorber has a la rge co l umn density (Nh Ant ~ 10^^ cm~^, 
iWozniak et"al] Il998t ISambruna et al.l ll999j) but unlike 
3C 111, the absorber in 3C 44 5 covers the source almost 
completely (/c ~ 0.8: ISambrun a et al. 1999). The reflec- 
tion fraction in 3C 445 is also very small < 0.2, 
but there is a very strong Fe Kg emission line with 
EW ~ 150 eV (|Wozniak et al.lll998f) : these authors pos- 
tulate that the Fe Ka emission line arises in a shell of cold 
material with Nh ~ (2-5) x lO^'^ cm~^ which is isotrop- 
ically irradiated from the central source. Thus, 3C 445 
and 3C 111 might in fact be quite similar, with the only 
difference being the covering factor of the absorber. 

The partial covering absorber model for 3C 111 would 
be best tested by UV and soft X-ray observations, be- 
cause these absorbing clouds should also exhibit emis- 
sion lines and absorption edges in these regimes. With 
combined UV and X-ray information, it would be pos- 
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sible to place constraints on the physical conditions and 
possibly the location of the absorbing material. Unfortu- 
nately, the Galactic column density is so large {NH,Gai ~ 
10^^ cm~^) that these observations are not possible. 
However, high S/N UV and X-ray observations of other 
BLRGs, which do not suffer from such a high Galactic 
extinction, might reveal that intrinsic absorption is not 
uncommon in BLRGs. If the covering factor of the ab- 
sorber is small, as inferred for 3C 111, then presence of a 
partial covering absorber could easily be overlooked be- 
cause the resulting continuum curvature is in the same 
region as the expected Fe Ka emission line. 

Alternatively, the partial covering model presented 
here might simply be a parameterization of the trun- 
cated accretion disk (refsch) model presented in i)5.3.1l 
in which significant continuum curvature near the 
Fe Kaline is present (Fig. (HJ. When a Gaussian line 
was added to that model in place of a disk emission line, 
the EW and a are very similar to those obtained when 
using the partial covering continuum (see Fig. Ill|l . Fur- 
thermore, when Seyfert Is that possess broad Fe Ka line 
features are fitted with a partial covering model , one ob- 
tains 7Vh,2 ~ 3-4 xl023cm-2 and fc < 0.35 llGelbordI 
1200,*^ . As with 3C 111, even when the data from 3-8 keV 
are excluded from the fit, the values of A'h,2 and fc are 
similar (Gelbord, private communication). These values 
of iVH,2 are somewhat larger than found for 3C 111, but 
as discussed in ^5.3.11 the curvature in the REFSCH model 
appears to be more extreme for larger values of r2/27r; 
when a simulated data set based on a refsch model with 
r2/27r = 0.75 and i?i„ — 6rg (using the fakeit command 
in XSPEC) is fitted with a partial covering model, we 
find A^H,2 ~ 3 X lO^'^cm"^. T hus it is possibl e that, as 
with the Seyfert Is studied bv lGelbordI l|2003D . the par- 
tial covering model for 3C 111 (Model #9) is simply a 
parameterization of continuum curvature and is not due 
to a real absorber. 

7. CONCLUSIONS 

In this paper, we present the results of an analysis 
of simultaneous observations of the BLRG 3C 111 with 
XMM-Newton and RXTE. The flux is moderately vari- 
able, but there is no evidence for significant spectral vari- 
ability. The continuum emission has been fitted with a 
wide variety of models and in all instances there are unac- 
ceptably large residuals at low energies. These are well 
fitted by a Gaussian component, which is included in 



all of the continuum models and is attributed to uncer- 
tainties in the calibration of the MOS detectors. Clear, 
broad, residuals are also present near the Fe Ka emission 
line region for all models. These can be fitted with a 
broad Fe Ka emission line with an equivalent width of 
~ 40-100 eV, which is weaker than those observed in 
Seyferts Is. The high energy RXTE data strongly disfa- 
vor the presence of a strong Compton reflection bump, 
even if the disk is highly ionized. This result gives strong 
support to the hypothesis that the geometry of the ac- 
cretion flow in BLRGs is different from that in Seyfert 1 
galaxies. 

We find that continuum curvature is a primary source 
of the broad residuals seen in the Fe Ka line region, al- 
though line emission is also required. The data are con- 
sistent with a model in which a weak Compton reflection 
bump and Fe Ka disk emission line are formed in a trun- 
cated accretion disk which transitions to a vertically ex- 
tended structure (such as an ion torus) at small radii. A 
less complex model, in which the primary X-ray source is 
partially covered by a dense (10^'^ cm~^) absorber, also 
provides a satisfactory fit to the data. In both models, 
an Fe Ka emission line that most likely arises in a distant 
reprocessor, such as a Compton-thin obscuring torus or 
dense clouds along the line of sight, is also present. Given 
these data, it is not possible to distinguish between these 
two models. However, it is very likely that the partial 
coverage model is simply a parameterization of the more 
complicated model in which the primary emission is re- 
processed by an truncated accretion disk. Moreover, the 
partial coverage model is unappealing in view of the small 
inclination angle of the jet in 3C 111. 
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APPENDIX 
INCLINATION ANGLE 

Following the methods of lEracleous et aP l)1996|) . we use the observed superluminal motion in the radio jet of 3C 111 
and the projected linear size of the radio lobes to place constraints upon the inclination angle of the jet, and thus the 
accretion disk. 

The measured proper motion, fi, in the jet is 1.54 ± 0.2 arcsec yr~^ ijVermeulen fc Cohenlll99'^ . The apparent 
velocity, relative to the speed of hght, is given by Papp — 47.4 fi z h~^, where z is the redshift, and h is Hubble's 
constant, in units of 100 km Mpc~^. Thus, fo r 3C 111 ffarm = 5.1 ± 0.7, assuming h = 0.7 , whi ch implies either 
i < 1 3° or 10° < i < 26°. We note that Er acleous fc Halpernl l|1998D . lR,evnolds et^ l|1998D and lEracleous et alJ 
1)20001^ neglected the factor of h in the denominator, and thus underestimated Papp and overestimated the upper limit 
on the inclination angle i. 

A lower limit on the inclination can be inferred from the size of the radio lobes. iNil^^yn ^1^ ^l.| il993 [) m easure the 
largest angular size of the radio lobes to be 275", which translates to 184.8 kpc/h. iGavazzi fc Pcrola (1978) find that 
for powerful radio sources (log P1215 mhz > 25.0), the radio lobes have an average intrinsic size of 180 kpc h^^, with a 
tail of sources which extends to 500 kpc h^^. Assuming the radio lobes of 3C 111 have an intrinsic size of 500 kpc h~^, 
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the inclination angle must be greater than 21.7°. However, it is possible that 3C 111 is a giant radio galaxy (GRG) 
and the intrinsic size of its radio lobes could be much larger. There are some enormous giant radio galaxies, such 
as 3C 236 (3 Mpc /i"^; Nilsson et al. 1993) and NVSS 2146+82 (2 Mpc ft.-^ Palma et al. 200 0), but these sources 
are r are and the size distribution of GRGs drops off rapidly for sizes greater than 1 Mpc l|Schoenmakers et alJ 
HqoT). Assuming 3C 111 is not larger than 1 Mpc h'^, we have i > 10.6°. Thus we adopt a conservative range in the 
inclination angle, 10° < i < 26°, keeping in mind that unless 3C 111 is a GRG, it is very likely that 21° < i < 26°. 
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Table 1. Observation Log 



Instrument Exposure Time (ks) Exposure Time (ks) Count Rate (counts s '^)^ 

post-filtering post-filtering 

XMM- Observation Date: 2001 Mar 14 13:03 - 2001 Mar 15 01:16 

pn 42.3 23.4 10.89 

MOS 1 44.0 27.8 3.22 

MOS 2 44.0 28.4 3.22 

RGS 1 44.6 43.5 0.08, 0.05 

RGS 2 44.6 42.5 0.10, 0.04 

RXTE - Observation Date: 2001 Mar 14 08:49 - 2001 Mar 17 03:36 

PCA 56.6 6.57 

HEXTE 18.2 1.14 

HEXTE 1 17.8 0.86 



''The count rate was measured for the same interval used for fitting. (For the pn data, 0.5-10 keV was used.) 
''The 1st and the 2nd order count rates, respectively. 



Table 2. Best Fit Continuum Parameters 



Model 



Model Number 



Model Parameters 



xVd.o.f. 



Powerlaw 



A'h= (8.1 ±0.1) X 10^1 cm-2 1005/895 
r = 1.72 ±0.01 



Broken Powerlaw 



Nh= (8.7 1[]-^) xlO^i cm- 
i 1 - 1.8b _Q 03 
£^broak = 2.9 ± 0.2 keV 
r2 = 1.69 ±0.01 



928/893 



Powerlaw + Compton 
Reflection (neutral) 
^agdzia^^^dziareki^95) 



Nh= (8.1 ±0.1) xlO^i cm 
P _ 1 74 +0.01 

' — J^- -0.02 
£^fold > 440 keV 



. n -|o +0.06 
"■^•^ -0.07 



999/892 



Powerlaw + Compton 
Reflection (ionized) 
(Magdziarz fc Zdziarski 1995: Done et al. 1992) 



N„= (8.1 ±0.01)xl02i cm" 
P _ 1 74 +0.02 
^ — J^- -0.01 
Efoid > 440 keV 
Q/27r = 0.13 ± 0.06 
^ < 20 erg cm s~^ 



997/890 



Powerlaw + Bremsstrahlung 



Nh 



:H^)xl02i a 

1.68 IH? 
kT = 0.9 ± 0.2 keV 



(8.5 : 
r = 



938/892 



Powerlaw + Gaussian 



= (8.0to2)xlO^^ 

r 

E -- 



0.21 

-- 1.69 ±0.01 
1.5 tol keV 
0.5±0'.l kcV 



916/892 



''The error bars are the 90% confidence interval for 1 degree of freedom (i.e. Ax^ = 2.706). 

'"The observed 2-10 keV flux, as measured with the MOS 1 and 2 cameras is (5.8 ± 0.1) X 10~^^ 
at the distance of 3C 111, is 3 X 10*"* erg s~^. 

'^The Bremsstrahlung component has a 0.5-2 keV unabsorbed flux of (5 ± 2)xl0~^^ erg s~^ cn 
0.5-2 keV band and 1% of the flux in the 2-10 keV band. 

"^The unabsorbed flux of the emission line 1.4 4 X 10"^'^ erg cm~^ s^^. 



^cm ^. The unabsorbed luminosity. 



and contributes 10% of the flux in the 
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Table 3. 
Models 



Best Fit Parameters for Combined Continuum and Fe Ka Line 



Model 



Model # 



Model Parameters ^ 



xVd-o.f. 



Powerlaw + Compton 
Reflection (refsch) 
+ Diskline 

Powerlaw + Compton 
Reflection (refsch) 
+ Diskline 
+ Narrow Gaussian 

Powerlaw + Compton 
Reflection (rf-pexrlv) 

Powerlaw + Compton 
Reflection (rf-pexriv) 
+ Narrow Gaussian 

Blurred Powerlaw + Compton 
Reflection (rf-pexriv)^ 
+ Narrow Gaussian 

Powerlaw 
+ Partial Covering 
+ Gaussian 



7a (see i|S.::i.lt 
7b (see 

8a (see A^.'i.'ll 
8b (see 

8c (see ^5:3:21 
9 (see ^5:3:31 



N„= (7.9 ±0.01) X 10- 
Bfoid < 120 keV n/2iT = 0.20 



Rin = 110 



+350^ 



cm ; 
+0.04 
-0.06' 



r = 1.66 ±0.02 

^ < 4 X 10^ erg cm s" 



;.4±0.1 keV; EWa 



32 



+6 



eV 



Nh= (7.8 t^;?)x 10^1 cm-^; F 



1.63 



+0.03 
-0.02 

Sfoid < 150 kcV; Q/2n = 0.09 to.ol; C < 5 x 10^ erg cm s 
Rin < 500 rg; Ea = 5.8 t^f, EWa 
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+6 



EWn 



25 ; 



+0.01 



Nh= (8.3±0.1)xl0^i cm-2; r = 1.64 
n/27r = 0.3 ± 0.1; ^ = 4 x 10^ erg cm s" 



Nh= (8.2±0.1)x102i cm-2; F = 1.66 
n/27r = 0.17 ± 0.06;a g = 4 I2 x 10^ erg cm s" 



+0.01 
-0.02 



EW„ = 15 



+12 



eV 



Nh-- 
n/27r = 



8.2 X 10-^^ cm 



1.65 



0.2; 5 = 3.5 X 10^ erg cm : 



EWn = 15 eV 



JVh= (8.2 tO-i)x 10^1 
A'h,2 = (1.5 j;?;I)xlO' 
: 6.4 ±0.2 keV; a = 0.2 



23 



-2; F 
-2 



+0.2 



1.76 
= 0.15 



+0.02 

-0.04 
+0.07 
-0.05 



keV; EW = 30 ± 20 eV 



1039/994 
1031/993 

1050/997 
1043/996 

1032/988 
1000/994 



''The soft Gaussian component (see i|5.1l Model #6) is included in all of these models. The flt parameters were allowed to vary, but the 
best-flt parameters are consistent with those found for Model #6: E ~ 1.5 keV; a = 0.5 keV; and an unabsorbed flux of ~ lO^^^erg s~^ cm'^. 
''The error bars correspond to the 90% confidence interval for 1 degree of freedom (i.e. Ax^ = 2.706). 

'^The diskline model has /3 = -3, i = 26°, and Rout = 500rg. The normalization of the emission line has been tied to Q/2n, such that 
EW = 160 X Q/2tt eV. The energy and equivalent width of the diskline are denoted by E^ and EWd- 
"^The narrow Gaussian has a fixed energy of 6.4 keV and width of 0.02 keV. The equivalent width of the narrow line is denoted by EWn- 
'The folding energy was restricted to be greater than 100 keV. 



SThe model was convolved with a diskline model with /3 = -3, i = 26°, Rin = 500rg and Ro 



; 400r„ 



RXTE PCA 2-20 keV 



- 


T 











5.0x104 1.0x105 1.5x105 2.0x105 2.5x10^ 

Time (s) 



Fig. 1. — The observed RXTE PCA lightcurve in the 2-20 keV energy band. The dashed vertical Unes indicate the interval during which 
the XMM data were taJjen while the dashed horizontal line indicates the average count rate over the entire observation. 
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Fig. 2. — Top panel - Observed 0.4—100 keV spectrum with the different continuum models described in i|5.1l overlaid (the fits excluded 
the data from 4.5-7.5 keV): Powerlaw (top left), Powerlaw + Compton Reflection (top right), Powerlaw + Bremsstrahlung (bottom left), 
and Powerlaw + Soft Gaussian (bottom right) models The RGS data have been excluded and the MOS have been binned for clarity. Lower 
panels - Ratio of the data to the model for, top to bottom, MOS 1, MOS 2, PCA, and HEXTE (cluster - crosses, cluster 1 - flUed circles). 
The position of 6.4 keV Fe Ko emission line at the redshift of 3C 111 is indicated by the vertical dashed line. The line is clearly seen in 
the PCA data, but is less obvious in the MOS 1 and MOS 2 data. 
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Flc;. 3. — Confidence contours (90%) in the Nh - T plane for four different continuum models: Powerlaw (solid); Powerlaw + Compton 
Reflection (dotted); Powerlaw + Bremsstrahlung (dashed): and Powerlaw + Soft Gaussian (dash-dot). The best-fit powerlaw index for the 
Powerlaw -|- Bremsstrahlung and Powerlaw -|- Soft Gaussian models is significantly steeper than that found for the Powerlaw and Powerlaw 
-|- Compton Reflection models. 
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Fig. 4. — Confidence contours in the Efoid - n/27r plane for the Powerlaw + Compton Reflection model. The inclination was fixed to 
26°. The reflection fraction is small, but non-zero. 
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Fig. 5. — Confidence contours (90%) for the Gaussian fit to the Fe Ka emission fine to the residuals from the Powerlaw + soft Gaussian 
continuum fit (model #6, solid). For comparison, we show the contours for the Gaussian fit to the residuals from the Powerlaw + Co mpton 
Reflection fit (model #3, dashed) as well as the simultaneous fit to Fe Kaand the continuum over the 3-10 keV interval (dotted, see il5.2i . 
The equivalent width is defined to be the line flux divided by the speciflc continuum flux at the rest energy of the emission line, E. 
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Fig. 6. — Confidence contours (90%) for the diskline fit to the Fe Ka emission fine in Rin—E space, using the Powerlaw + soft Gaussian 
continuum model (Model #6). The shape of the line profile is sensitive to the inclination i and the powerlaw emissivity index, P, thus we 
present a set of Ri„-E contours for i = 26° (solid), 18° (dashed), and 10° (dotted) using two different values of /3. 
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Fig. 7. — Confidence contours (90%) for the diskUnc fit to the Fe Ko emission Une in Rin—E space, using the Powerlaw + Compton 
Reflection continuum model (Model #3) for i = 26° (solid), 18° (dashed), and 10° (dotted) and /3 = -2.5. 
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Fig. 8. — Ratio of a simulated PCA data set based upon model #7a (without the disk emission line component) to the Powerlaw + 
soft Gaussian model (Model #6). The shape of the residuals is indicative of curvature in the former continuum model, which can mimic a 
broad Fe Ka line. 
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Fig. 10. — Confide nce co ntours in tlie column density and covering fraction for the partial covering powerlaw model described in 
t|5.3.3l As described in i|5.3.3l the absorber is expected to produce an Fe Ko emission line. Overlaid on the contours are lines of constant 
Fe Ko equivalent width (10,15, and 20 eV) as a function of A'^h,2 a^nd fc as prescribed by equation (1). 
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Fig. 11. — Confidence contou rs for the Gaussian fit to the residual Fe Kq emission fine from the partial covering model (with the soft 
Gaussian included) described in H5.3.3I The equivalent width is defined to be the line flux divided by the specific continuum flux at energy 
E. 



